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Abstract 


This paper presents a study on the effect of inhomogeneous compression of gas diffusion layer (GDL) caused by the channel/rib structure of 
flow-field plate. The experimentally evaluated properties are GDL intrusion into the channel, gas permeability, in-plane and through-plane bulk 
electric conductivities, and contact resistances at interfaces as a function of compressed thickness of GDL. It was found that the GDL is compressed 
very little under the channel whereas GDL under the rib is compressed to gasket thickness. The compression of GDL reduces gas permeability and 
contact resistance, and improves bulk conductivity. Hence, the inhomogeneous compression of GDL may lead into significant local variation of 
mass and charge transport properties in the GDL. These effects have been ignored in most of the published modeling studies. This contribution, 
part I, covers the experimental setup and measurement results, and a model which takes the inhomogeneous compression of GDL into account is 


presented in part II. 
© 2006 Elsevier B.V. All rights reserved. 
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1. Introduction 


Proton exchange membrane (PEM) fuel cells are electro- 
chemical devices that convert the chemical energy of reactants 
directly into electrical energy. This technology enables high effi- 
ciency and energy density compared to conventional internal 
combustion engines, thereby making the technology attractive 
for automotive, portable, and stationary applications. Further- 
more, the only exhaust from PEM fuel cells is water, which 
makes them favorable from the environmental point of view. 

PEM fuel cell consists of bipolar plates with channels 
machined on either side for reactant distribution over the elec- 
trode surface, membrane electrode assembly (MEA) where the 
electro-chemical reactions and proton transport takes place, 
and porous gas diffusion layers (GDL) sandwiching the MEA. 
Among the components of PEM fuel cell, much research effort 
has been put on catalysts, membranes and bipolar plates. 
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Until recently, less attention has been paid to GDL even 
though it plays an important role in fuel cell operation. The 
main functions of GDL are to provide a passage for reactant 
access and product water removal, to conduct electricity and heat 
between adjacent components, and to provide mechanical sup- 
port for the MEA. These functions impose requirements on the 
electrical and mechanical properties of GDLs, i.e. high gas per- 
meability and suitable water management properties, electrical 
and thermal conductivity, and chemical and physical durability. 
Typically, GDLs are made of highly porous carbon-fiber based 
paper or cloth. 

The properties of GDL and interfaces between the GDL and 
the electrode and flow field plate are strongly dependent on com- 
pression pressure. Generally speaking, increasing compression 
improves electrical and thermal conductivity, and impedes reac- 
tant transport and water removal. A flowchart of the effects of 
increasing compression is presented in Fig. 1. 

Pressure onto the GDL is exerted by the flow field plates, 
which usually feature flow channels for reactant distribution. 
The rib/channel-structure creates an inhomogeneous compres- 
sion distribution, because the GDL under ribs is compressed to 
gasket thickness while GDL under the channels remains mostly 
uncompressed and intrudes into the channel, see Fig. 2. 
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Fig. 1. The effects of increasing compression of GDL. 


Inhomogeneous compression distribution leads into spatial 
variation of GDL and interface properties. Fluid permeability is 
lower in the GDL under the ribs due to loss of porosity, which 
increases mass transport overpotential in those areas compared 
to areas under the channel. Similarly, electrical and thermal bulk 
conductivities are improved and contact resistances at the inter- 
faces between the GDL and electrode and flow field plate are 
smaller than those under the channel. 

These effects exist in all fuel cells with normal flow field 
plates but cannot be measured directly due to small scale of the 
phenomena. Thus, the only available option is to characterize 
the GDLs ex situ and model the effect. Furthermore, most of 
the published fuel cell models do not account for the inhomo- 
geneous compression of the GDL and its effects. Usually GDL 
thickness, porosity, contact resistances and conductivities are 
assumed constant over the cell area. There may be a significant 
discrepancy between the modeled results and practical situation 
due to these assumptions. Only few studies, such as by Sun et al. 
[1] and Zhou et al. [2], were found by authors, in which the effect 
of inhomogeneous compression of GDL is taken into account. 
However, Sun et al. did not account for the contact resistance, 
and Zhou et al. ignored the effect of compression on the bulk 
conductivity of GDL. 


Bipolar plate 
GDL intruded into channel 


GDL compressed under rib 


Fig. 2. Schematic illustration of GDL deformation under compression. The 
GDL partially intrudes into the channels and the parts under the ribs are more 
compressed. 


Effects of compression and GDL properties on fuel cell per- 
formance have been studied by several groups, e.g. Lee et al. 
[3], Ge et al. [4], and Ihonen et al. [5]. The results show that the 
compression force and physical properties of GDL must be con- 
sidered together and there is an optimum compression pressure 
and compressed thickness for each GDL which trades off the 
competing demands of mass, charge and heat transport within 
the GDL. 

GDL parameters, such as permeability, both electric and ther- 
mal conductivity and contact resistance between components 
have been investigated, e.g. Dohle et al. [6], Williams et al. [7], 
Feser et al. [8], Ihonen et al. [9], Mishra et al. [10], Wang and 
Turner [11], Vie and Kjelstrup [12] and Khandelwal and Mench 
[13]. However, in contact resistance studies [9—11], the effect 
of compression on bulk conductivity was ignored. Furthermore, 
Thonen et al. [9] and Natarajan and Nguyen [14] found that fuel 
cell operating parameters or experimental operating conditions 
also affected the contact resistance. Relating to compression, it 
has been observed that excessive compression damages the car- 
bon fibers in GDL materials, e.g. Wilde et al. [15], Escribano et 
al. [16] and Matsuura et al. [17]. 

Some of the actual GDL properties have been included in 
published fuel cell models. Chu et al. [18] considered the non- 
uniform porosity due to presence of water in GDL. Berning and 
Djilali [19], Natarajan and Nguyen [20], Inoue et al. [21], and 
Jang et al. [22] modeled the mass and charge transport with var- 
ious operating parameters and studied the impact of geometric 
parameters of flow field and material properties of GDL such 
as porosity and thickness. Pharoah [23] applied anisotropic gas 
permeability into their models, and Meng and Wang [24] and 
Pharoah et al. [25] modeled the effect of anisotropic electrical 
conductivity and observed large variations between isotropic and 
anisotropic cases. However, the authors did not find any models 
that took into account the spatial variations in permeability, and 
electrical and thermal conductivity due to inhomogeneous com- 
pression. This may cause significant errors in modeled results, 
e.g. in the prediction of current density distribution, and there- 
fore it is worth studying how the compression pressure affects 
charge and species transport. 
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The purpose of this contribution is to experimentally eval- 
uate the GDL parameters (GDL intrusion into channel, gas 
permeability, electrical bulk conductivities and contact resis- 
tances) as a function of compressed thickness. Determination 
of thermal conductivity and thermal contact resistance was left 
outside this work. The evaluated values will be used in a model- 
ing study, in which the effects of inhomogeneous compression 
on the local phenomena are taken into account. These results 
of modeling studies are given in part II of this contribution 
[26]. 


2. Experimental 


A commercial gas diffusion media, SGL SIGRACET 10- 
BA carbon paper (made by SGL Carbon Group), was used in 
the measurements. This media has uncompressed thickness of 
380 um, is treated with 5 wt%. PTFE for wet proofing, has no 
microporous layer, and has relatively high open porosity of 88% 
[27]. 


2.1. GDL intrusion into channel 


The amount of GDL intrusion and shape change of GDL 
under compression must be known in order to evaluate the 
local parameter values of GDL. The initial GDL thickness and 
gasket thickness determine the deformation of GDL regardless 
of the compression pressure when gaskets having low com- 
pressibility are used. Hence, the measurements are conducted 
from the perspective of thickness, not of compression pres- 
sure. 

The experimental setup is illustrated in Fig. 3. In order to con- 
trol the thickness of compressed GDL, steel gages were inserted 
between two aluminum plates. The thickness of steel gages was 
varied between 150 and 350 pm. Steel gages were placed near 
the GDL, having an area of 1.6 cm?, to allow precise control of 
the compressed GDL thickness. A channel with a width ranging 
from 0.6 to 2 mm was machined into the bottom aluminum plate. 
The parts were clamped together with four 8 mm bolts, using a 
torque from 0.8 to 1.5 N m, depending on the desired thickness. 
GDL intrusion was measured with a dial indicator (ND 221B by 
Heidenhain Corporation). 
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Fig. 3. Schematic illustration of GDL intrusion measurement setup. 
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Fig. 4. Gas permeability measurement setup: (a) top view, (b) cross-section 
view, and (c) enlarged cross-section view. 


2.2. In-plane gas permeability 


The method of gas permeability measurement in this paper is 
fundamentally similar that used by Mikkola [28], but modified 
to allow the control of the thickness of GDL. The top view of 
the measurement device is illustrated in Fig. 4(a), and its cross- 
section is illustrated in Fig. 4(b) and (c). A circular GDL piece of 
30 mm diameter was mounted concentrically on top of the inlet 
in the steel base plate. Four steel gages were situated in each 
corner of the base plate in order to enable the precise control 
of the GDL thickness. The upper base plate presses and fastens 
the inserted GDL and steel gages with eight bolts to ensure the 
assembly. The radii of air inlets, r1, and of GDL, r2, were 8 and 
15 mm, respectively. Gas pressure at inlet rı was pı and the gas 
pressure p2 at r2 was ambient pressure, since air was discharged 
freely into the surroundings. 
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The measurements were conducted by varying the com- 
pressed thickness of GDL from 150 to 350 wm with 50 pm steps 
controlled by steel gages. Dry air was used as the flowing fluid. 
The flow rate of air was varied in the range of 10-310 cm? min7! 
controlled by a mass flow controller (Model 5850S by Brooks 
Instruments). 

The porosity loss in GDL under compression increases the 
pressure drop and reduces mass transfer. The flow resistance is 
characterized by the gas permeability, which can be estimated 
by measuring pressure difference between inlet and atmosphere 
and by applying Darcy’s law: 


k 
v=—-Vp (1) 
n 
where vis the flow velocity, k the gas permeability, 7 the dynamic 
viscosity and p is the pressure. In cylindrical measurement 
geometry the flow can be assumed to be radially symmetric. 
By applying conservation of mass and ideal gas law, one can 
deduce the permeability from Eq. (1): 


k= RT n Vi l r2 
~ Voh age r 
Pi — P? 1 


RT n Vi In F3 (2) 
~ aVoh 2_ p2 r 
0 A (pm + P2) — p5 1 


where V; is the source volume flow rate, R the gas constant, T 
the temperature, h the cylinder height (here equivalent to com- 
pressed GDL thickness), and Vo is the standard molar volume. 
The pressure difference was measured with a manometer (Type 
MM3K by HK Instruments Oy). 


2.3. Through-plane GDL conductivity 


The conventional experimental setup to evaluate the contact 
resistance between GDL and current collector typically consists 
of two graphite (or steel) current collectors, two supporting base 
plates, sandwiching one GDL as illustrated in Fig. 5(a), see, e.g. 
[10,29]. 

In the previous studies, various compression pressures were 
exerted on both sides of the base plates, which led into a change 
in the thickness of GDL and thus also in bulk and contact resis- 
tances. With these conditions, total resistance of the system was 
measured using four-point probe method. The measured total 
resistance Rzmeas,conv(z) with this conventional experimental 
setup is expressed as 


Rz, meas,conv(Z) = 2RGR + Rb,Gpi(Z) + 2Rc,GR/GDL(Z) (3) 


where Reor and Rp,Gpi(z) denote the bulk resistances of graphite 
current collector and GDL in through-plane direction, respec- 
tively. Re,GR/GDL(Z) is the contact resistance between graphite 
current collector and GDL. 

The ReGr/GpL(z) is typically evaluated by neglecting RaR 
and with a simple subtraction of Rp,GpL(z) which is based on 
the value given by manufacturer. In this conventional method, 
Rb, GDL (Z) is considered constant and the effect of compression 
on the bulk conductivity of GDL is simply disregarded. In fact, 
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Steel gage 
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Zz 
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Fig. 5. Contact/bulk resistance measurement setup: (a) conventionally applied, 
and (b) system used in this study. 


the bulk conductivities of GDL strongly depend on the compres- 
sion pressure as discussed in Section 3.3, and consequently this 
conventional method inherently contains errors in the evaluated 
contact resistance values. 

The compression pressure affects the bulk resistance of GDL 
and can be expressed as a function of compressed GDL thickness 
g: 

Roc.) = — (4) 
iii 0z, GDL(Z)A 

where oz pL is the through-plane conductivity of GDL which 

is, in fact, a function of z, and A is the area of GDL (1.6 cm? in 

the measurements). 

As shown by Nakamura et al. [29], it is intrinsically difficult 
to evaluate the bulk GDL resistance separately from measured 
total resistance since two variables (Re,GR/GDL(2), Rb,GDL(Z)) 
are included in one Eq. (3), both of which are functions of the 
compressed thickness. Cunningham et al. [30] successfully iso- 
lated the bulk resistance of bipolar plate by their calibration 
procedures, but they could not separate the bulk and contact 
resistances of carbon paper. 

Here, an effort is put to evaluate Ry GpL(z) separately from 
the measured resistance, which allows investigation of the com- 
pression effect on GDL bulk conductivity alone. In this study, 
two to five pieces of GDLs were placed on top of each other 
as illustrated in Fig. 5(b). Non-conductive tapes were pasted on 
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current collectors in order to electrically isolate the end plates 
form each other. The thickness of each GDL was controlled by 
steel gages having thickness from 150 to 350 wm. For example, 
when three GDLs were used and each thickness of GDL was 
assumed to be 250 um, the sum of thicknesses of steel gages 
and non-conductive tapes was set to 750 wm. This means that 
no matter how many GDLs were used, both Re,GR/GDL(Z) and 
Rp,GpL(z) were unchanged when the compressed thickness of 
GDLs was set constant. In this experimental setup, the total 
resistance Rz meas(n) which is a function of number of GDLs, n, 
can be expressed as 


Rz meas(Z, n) = 2RoR + n Rb GDLC) + 2Re,GR/GDL{) 
+(n — 1)Re,GDL/GDL(Z) (5) 


where Re GpL/GpL is the contact resistance between two facing 
GDLs. If this Re GDL/GDL can be eliminated from Eq. (5), it 
reduces into 


Rz,meas(Z, n) = 2Rar + Rp,cpi(Z) + 2Re,GR/GDL(Z) (6) 


Thus, when one plots Rz meas(n) as a function of number 
of GDLs, the slope of the graph indicates the Rp,GpL(z). The 
detailed explanation about the derivation of the Re GryGpL(Z) is 
explained in Section 2.6. 

In order to eliminate Re GDL/GDL, highly conductive silver 
particles were sputtered onto GDLs. This process consisted of 
vacuum evaporation of coating material and deposition onto the 
target GDL. The deposition time was chosen so that the thickness 
of the silver coatings corresponds to 150 nm layer. This method 
was previously shown to successfully reduce the contact resis- 
tance [31,32]. According to Caillard et al. [33], the diffusion 
depth of deposited particles is in the range below 100 nm, which 
is small enough compared to the thickness of GDL in order 
not to affect the bulk conductivity significantly. This allowed 
assuming that the contact resistance between facing GDLs 
was negligibly small without changing the bulk properties of 
GDL. 

All GDL surfaces that came into contact with other GDLs 
were silver coated, but GDL surfaces facing the graphite 
current collectors were left untreated. Therefore, contact resis- 
tance between the GDLs could be assumed negligible and 
Re,GR/GDL(Z) could be evaluated. Resistance measurements were 
conducted by applying a current through the sample with 50 mA 
steps from 50 mA to 1 A, and simultaneously recording the volt- 
age drop in order to calculate the total resistance. 


2.4. In-plane GDL conductivity 


The experimental setup used for investigation of in-plane 
conductivity as a function of compressed thickness of GDL is 
illustrated in Fig. 6. Fig. 6(a) illustrates the top view of the setup. 
GDL is placed on the base plate and both ends are compressed 
by graphite current collectors clenched with two 8mm bolts. 
Steel gages were inserted between graphite current collectors 
and base plate, and thus the GDL thickness under the graphite 
current collector was always fixed to the steel gage thickness 
of 250 wm. The separation between two current collectors was 


Voltage measurement 


Base plate 


) 


Zz 
(b) n 


Fig. 6. In-plane conductivity measurement setup: (a) top view without plastic 
plate, and (b) cross section view. 


varied in the range from 1 to 32mm. Fig. 6(b) illustrates the 
cross section view of the setup. A non-conductive plastic plate 
was placed between the graphite current collectors. This plate 
allowed the accurate control of the compressed thickness of GDL 
underneath it and maintained the separation. The plastic plate 
was clamped to the base plate with two 8mm bolts, and the 
thickness of GDL under plastic plate was controlled from 150 
to 350 um with steel gages. 

A current in the range from 50 mA to 1 A with 50mA step 
was applied and the voltage drop was measured to calculate the 
total resistance. The measured total resistance Ry meas(x,z) con- 
sists of the bulk resistance of graphite current collector RoR, 
bulk resistance of GDL under the graphite current collector 
Rb ,cpL (250) and under the plastic plate Rp.GpL(%z) and the 
contact resistances between graphite current collector and GDL 
Re,GR/GDL* 


Rx, meas(%, Z) = 2RGR + Rv,Gpi(250) + Rv,Goi@, z) 
+ 2Re,GR/GDL (7) 


In order to separate the bulk GDL resistance from the total 
resistance and to investigate the effect of compression, follow- 
ing facts and assumptions should be underlined. First, since the 
steel gages under graphite current collector fixed the compressed 
thickness of GDL to 250 um, both ReGricpL and Rp,cpL(250) 
remained constant through the measurement series. Second, 
the conductivity of graphite current collector was at least one 
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order of magnitude larger than that of GDL, and therefore the 
current path within graphite current collector was assumed to 
concentrate to the edge of the plate and Ror remained con- 
stant. And third, since the separation between the graphite plates 
was at least an order of magnitude larger than the thickness 
of compressed GDL, it was assumed that the effect of current 
flowing in through-plane direction from the graphite plates was 
negligible, allowing us to study only the changes in in-plane 
direction. 

Based on the above facts and assumptions, Rp,GpL(,z) can be 
derived by subtracting Ry meas(x) from the other measured values 
with different separations between graphite plates. For example, 
with two different separation, xı and x2 (x; <x2), and with same 
thickness z1, the bulk GDL resistance can be calculated from the 
subtraction of total resistances as 


Rx, meas(X2; 21) — Rx,meas(41, 21) = Rb,GoL@2 — x1, 21) 
= Rb GDL(Ax, z1) (8) 
Then, the in-plane conductivity ox gpL(z) can be calculated 
from the bulk resistance as 
1 Ax 


= 9 
Ox GDL(Z) one a (9) 


where Ax, z and w are the current collector separation, com- 
pressed thickness and width (1 cm in all of the measurements) 
of the bulk GDL underneath the plastic plate, respectively. In 
order to increase the statistical accuracy, measurements were 
conducted using various current collector separations and the 
average value was used for the in-plane conductivities for each 
compressed thickness. 


2.5. Conductivity of other components 


The conductivities of graphite current collector and elec- 
trode on MEA were also evaluated by applying four-point 
probe method. This was done to calculate the contact resistance 
between GDL and electrode on MEA as explained in Section 2.4. 

For determining the conductivity of graphite, a graphite block 
having dimensions of 8mm x 14mm x 90 mm was used. A cur- 
rent lead was connected to both ends of the block and current was 
applied in the range from 50 mA to | A. The potential drop was 
measured at various positions of the block with distances from 
9 to 63 mm. The conductivity of graphite was solved from the 
subtraction of bulk resistances obtained with various lengths in 
similar fashion to the in-plane conductivity of GDL as described 
in Section 2.4. 

In the evaluation of the electrode conductivity, the measure- 
ment setup was fundamentally the same as illustrated in Fig. 6, 
but the GDL was replaced with a piece of MEA (PRIMEA 
series 58 by W.L. Gore & Associates, Inc.). Since the con- 
ductivity of electrode was expected to be smaller than that of 
GDL, the separation between two graphite current collectors 
was varied from 0.5 to 10 mm and the thickness of electrode on 
MEA was assumed to be 10 wm. Then the same procedure as 
for the in-plane conductivity of GDL was applied to evaluate the 
conductivity of the electrode. 


2.6. Contact resistance between graphite current collector 
and GDL 


Finite element method based simulation was employed to 
fit the contact resistance between graphite current collector and 
GDL using acommercially available PDE solver program COM- 
SOL Multiphysics 3.2a, since it was believed to be difficult 
to separate the contact resistance from total resistance exper- 
imentally. The conductivity of graphite was only an order of 
magnitude higher than that of GDL, and due to its larger thick- 
ness, the resistance of current collector also affected the results. 
Thus, positioning of current leads affected the profile of current 
going through the system, and it could not be assumed constant 
throughout the whole contact area of the components. There- 
fore, the actual measurement geometry was implemented as a 
3D modeling domain. Model domain was same as the experi- 
mental setup for through-plane GDL conductivity (see Fig. 5) 
which was composed of two graphite current collectors and GDL 
sandwiched between them. The GDL thickness was changed 
from 150 to 350 um in the modeling domain corresponding to 
the experimental study and the evaluated bulk conductivities of 
GDL and graphite were embedded. 

The electric potential in graphite, pgr, and GDL, gcgpL, were 
solved in the model. The governing equation for current density 
follows from the conservation of charge: 


V-~ = V(—0z,crV¢cr) = V(—9z,GpLV capi) = 0 (10) 


Measured potential differences were used as boundary con- 
ditions in the model where voltage leads were attached (see 
Fig. 5). 

Since there is a contact resistance between graphite current 
collector and GDL, the potential profile is discontinuous. This 
potential drop at the interface can be expressed with current 
density and the contact resistance. At the graphite current collec- 
tor/GDL interface, Neumann boundary condition was applied: 
-A(—oz, or Voor) = “CR — POPE a1) 

Re,GR/GDL 

where ñ is the normal vector of the interface. After solving the 
potentials, the total current, J, passing through the system was 
calculated by integrating the local current density. Re,GR/GDL 
was evaluated by changing its value in the model until the total 
current corresponded to the measurement result. The same mod- 
eling procedure was applied to a system that consisted of various 
numbers of GDLs in order to find an error estimate for contact 
resistance. 


2.7. Contact resistance between GDL and electrode on 
MEA 


An attempt was made to evaluate the contact resistance 
between GDL and electrode, Re GDL/ÆLE, With fundamentally 
similar measurement setup and procedure as introduced in Sec- 
tion 2.4. The difference was that the GDL was replaced with 
MEA and there was GDL pieces between the graphite current 
collectors and MEA. The separation between two graphite cur- 
rent collectors was changed in the range from 0.5 to 5mm. The 
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potential profile was modeled with similar domain of in-plane 
experimental setup (see Fig. 6) by applying the measured total 
resistance and conductivities of graphite current collector, GDL 
and electrode, and contact resistance between GDL and graphite 
current collector evaluated previously. Now the only unknown 
parameter was the contact resistance between GDL and elec- 
trode, and its value was varied in order to fit the modeled current 
with the measured value similarly as in Section 2.6. 


3. Results and discussion 
3.1. GDL intrusion into channel under clamping pressure 


The results from the measurements of GDL thickness under 
the channel, i.e. the sum of the measured GDL intrusion and 
steel gage thickness, are shown in Fig. 7. The initial thickness 
of GDL (380 um) is also included in Fig. 7 for comparison. 
It was expected that the channel width may cause a change in 
GDL intrusion since the mechanical force applied onto GDL is 
dependent on the geometry of adjacent components (channels 
and ribs). However, the results indicated that the thickness of 
GDL under the channel is not strongly dependent on the channel 
width and/or the thickness of the steel gages, and it remains 
almost uncompressed. For example, if a 250 wm thick gasket 
having low compressibility is used, the GDL is compressed at 
maximum 10 pum under the channel regardless of the channel 
width. 

Since the change in thickness of GDL caused by compression 
is mainly attributed to loss of pore volume, it naturally has an 
effect on mass and charge transport through the GDL. The fact 
that the GDL is practically not compressed under the channels 
leads into significant local variations in the gas permeability of 
the GDL, as well as electric bulk and contact resistances. This 
will be discussed in the following sub-sections. 


3.2. In-plane gas permeability 
Fig. 8 shows the calculated in-plane gas permeability of 


GDL as a function of compressed thickness. The gas perme- 
ability decreased non-linearly when the thickness of GDL was 
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Fig. 7. Thickness of GDL under channel. 
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Fig. 8. In-plane gas permeability as a function of compressed thickness of GDL. 


decreased by compression. The reduction in the permeability 
was as much as one order of magnitude when the GDL was 
compressed to approximately to 65% (250 um) of the initial 
thickness (380 um). The range of results is in good agreement 
with those published in literature for typical carbon paper type 
GDLs. For example, Mathias et al. [34] determined the in-plane 
permeability approximately to be in the range from 5 x 107!” to 
10 x 10~!* m? for Toray TGP-H-060 carbon paper compressed 
to 75% of initial thickness. 

In the range where the measurements were conducted, com- 
pression of the GDL mainly leads into loss of pore volume, 
not compression of bulk material. Therefore, porosity can be 
assumed to correlate directly with compressed GDL thickness. 
However, most of the previous modeling studies may not directly 
reflect the realistic situation occurring in a fuel cell since the two 
above-mentioned parameters were studied separately and their 
combined effect was basically not considered. In fact, the com- 
pressed thickness (and thus porosity) and gas permeability of 
GDL are closely linked to each other as shown in Fig. 8, and it 
is the combined effect that determines the mass transport within 
the GDL. 

Except for a study by Dohle et al. [6], none of experimen- 
tal investigations on the effect of compressed thickness of GDL 
on mass transport was conducted, whereas numbers of studies 
focused on the initial thickness, see e.g. [35-37]. In these stud- 
ies, the impact of the initial thickness on such parameters as 
sensitivity to water accumulation, and mass and charge trans- 
port resistance is fully addressed. However, these parameters 
are more closely associated with the compressed thickness of 
the GDL rather than the initial thickness. While the estimation 
of gas permeability from the porosity loss as the ratio of com- 
pressed and initial thickness is perceptive, one must note the fact 
that gas permeability is a non-linear function of the compressed 
thickness. 


3.3. Through-plane and in-plane GDL conductivity 


Fig. 9(a) and (b) illustrates the measured total resistance in 
through-plane measurement setup as a function of the number of 
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Fig. 9. Measured total resistance: (a) in through-plane direction as a function of 
number of GDLs, and (b) in in-plane direction as a function of length of GDL. 


GDLs and the measured total resistance in in-plane measurement 
setup as a function of GDL length, i.e. the separation between 
current collectors, respectively. The results from through-plane 
measurements are shown for GDL thicknesses 150—300 um for 
the clarity of the figure. Although the results for 350 um behaved 
similarly as the other curves for different thickness, the total 
resistances were significantly higher and thus were not incorpo- 
rated in figure. It is obvious from the results that the measured 
resistances are linear functions in both Fig. 9(a) and (b), and 
thus the in-plane and through-plane conductivity can be calcu- 
lated from the linear fits (see Sections 2.3 and 2.4). The resulting 
values for the in-plane and through-plane conductivity of GDL 
as a function of compressed thickness are illustrated in Fig. 10. 
The achieved values are averages from several repeated mea- 
surements (S—8 measurements for each point), and the errors 
correspond to the standard deviations. 

The conductivity in in-plane direction was three to four times 
higher than that in through-plane direction, which is not as large 
a difference as stated in the literature, see e.g. [34]. The pre- 
viously reported data for GDL conductivity varies a lot even 
among paper type GDLs. The reported values range from 5000 
to 23,000 S m7! [5,7,34,38] for in-plane conductivity and from 
300 to 1400S m7! [14,27,29,34] for through-plane conductivity. 
These variations can most probably be attributed to the physical 
and structural properties of different GDLs such as hydrophobic 
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Fig. 10. In-plane and through-plane conductivities as a function of compressed 
thickness of GDL. 


agent content, fiber diameter, density, areal density and imple- 
mentation of micro porous layer. Furthermore, the significant 
difference in the reported conductivities may arise from the 
drawbacks in measurement setups and inadequate consideration 
of contact resistance as discussed in Section 2.3. 

Both in-plane and through-plane conductivities monoton- 
ically increase as the compressed thickness of GDL was 
decreased. For example, when the GDL is compressed to a thick- 
ness of 250 wm, the in-plane conductivity increases 160% and 
through-plane conductivity over 10-fold compared to the values 
of uncompressed GDL. It is interesting to note that the con- 
ductivities have practically a linear dependence on the GDL 
compressed thickness. This is possibly because of the reduc- 
tion of porosity of GDL, which leads into shorter distances 
between the conductive carbon fibers and also into more and 
better contacts between the fibers. In an assembled cell, there 
are significant variations in the GDL thickness between the 
areas under the flow channels and under the ridges. Therefore, 
the assumption of constant conductivity may distort modeling 
results significantly. 


3.4. Contact resistances 


The area-specific contact resistance between graphite current 
collector and GDL simulated by applying the conductivities pre- 
sented in the previous sub-sections is given in Fig. 11. Error 
estimates were obtained by conducting the simulations again 
with considering variations in measured voltages and calculated 
conductivities. The contact resistance changed exponentially 
as a function of compressed thickness of GDL, as repeatedly 
reported in the literature, see e.g. [9,10,34]. Contact resistance 
values decreased down to two orders of magnitude from very lit- 
tle compression (350 um) to high compression (150-250 wm). 
This result is most probably due to the fact that the actual con- 
tact area of the GDL at the interface increased with compression 
pressure. 

There is a vast amount of literature discussing the con- 
tact resistance values, varying a lot in the range of 1-50 mQ 
cm? with moderate compression pressure (0.5-2 MPa) 
[2,9,10,14,29,34,38-41]. This large variation and difference 
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Fig. 11. Contact resistance as function of compressed thickness of GDL. 


from the values achieved in this study may be attributed to the 
fact that in the previous studies, the bulk resistance of the GDL 
and current collector was typically underestimated or simply dis- 
regarded, which in turn lead into overestimation of the contact 
resistance. The bulk resistance of current collector particularly 
has been assumed to be negligibly small because of its high 
conductivity in comparison to that of GDL materials, which 
also is often the case. For example, the conductivity of graphite 
used in this measurement was found to be 69700 + 300S m7! 
which was much higher than that of GDL (see Fig. 10). How- 
ever, the thickness of current collector is also almost one order 
of magnitude larger than that of GDL because of the structure of 
flow field, which is typically constrained to the millimeter scale 
[42,43]. Thus, one cannot doubtlessly assume that the voltage 
drop in graphite current collector is small enough to be neglected 
in these measurements. Hence, it was worth investigating how 
significantly the bulk resistance of graphite affected the charge 
transport. 

The electric potential profile achieved from the simulations 
is plotted from top to bottom at the center of x—y plane of the 
modeled domain (the through-plane measurement system, see 
Fig. 5) in longitudinal direction, as depicted in Fig. 12(a). In 
order to be able to distinguish the potential drop at the interfaces 
and in the bulk GDL, the enlarged potential profile in the region 
around the interface is illustrated in Fig. 12(b). According to 
Fig. 12(a), it is obvious that the gentle potential drop inside 
1 cm thick graphite current collector at both ends (2A@ggr) is 
comparable to that caused by contact and bulk GDL resistances. 
Itis worth noting that especially at the highest GDL compression 
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Fig. 12. Potential profiles within the model domain in longitudinal direction: 
(a) whole system, and (b) enlarged around the GDL. 


2Ag@GrR accounts for the dominant part of total potential drop. 
Hence, the resistance of current collector must be also taken 
into account properly when evaluating the contact resistance 
values in order not to distort the results. This is especially the 
case if the current collectors used in the measurements have 
lower conductivity than graphite, such as carbon compounds or 
polymer/graphite composites. 

The potential drops caused by different factors and their 
relative significance are summarized in Table 1. The ratios of 
potential drop caused by two contact resistances and bulk GDL 
resistance to measured total resistance A@meas, 2AGrR, / AYmeas 
and AgcpL/A@meas, respectively, are also tabulated in Table 1. 
Both 2Ag¢r, and 2A¢r./AGmeas increased as the compressed 


Table 1 
Predicted voltage drops at bulk and interfaces, and percentage of them from the total measured voltage drop 
Thickness (wm) Potential drop (mV) Ratio (%) 
2AGr, APGDL 2A¢cGR 2AGrR, /APmeas APGDL/ A Pmeas 
150 0.17 0.44 0.67 8.8 23.3 
200 0.42 0.74 0.67 17.1 30.4 
250 1.23 1.12 0.67 33.6 30.7 
300 2.35 2.24 0.67 39.4 37.6 
350 6.90 6.04 0.67 51.2 44.9 
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thickness of GDL was increased. However, it should be noted 
that 2A@r,/A@meas Was unexpectedly small, whereas many 
sources claim that the contact resistance accounts for the major 
part of the total resistance, see, e.g. [9,10,29]. The results pre- 
sented here highlight the fact that the contact resistances cause 
the major contribution to the total resistance only when the GDL 
is very little compressed. 

These results also indicate that the bulk resistance of the 
GDL can make a significant contribution into the total resis- 
tance, while contact resistance was found to be smaller than 
previously reported. Furthermore, the effect of compression on 
bulk resistance cannot be neglected in the evaluation of the con- 
tact resistances. Since the bulk conductivities of GDL changed 
linearly with compressed thickness as shown in Fig. 10, the 
method used in previous studies, in which bulk conductivities 
are assumed constant may not yield proper values for contact 
resistance. 

No literature data exists to verify the contact resis- 
tance between GDL and electrode on electrolyte membrane, 
Re GDLÆLE, Which has typically been simply ignored in mod- 
eling. However, it can be expected that Re GDLÆLE may be 
significantly higher than R..GryGpL since the electrode is com- 
posed also of non-electrically conductive polymer. Discussion 
by Ihonen [39] implies that the value may be roughly twice 
the value of contact resistance between GDL and graphite. 
Re,GDL/ELE İs an important parameter since a high contact resis- 
tance may consequently lead into hot-spots that enhance the 
membrane degradation and can severely decrease the life time 
of fuel cells. 

The Re GpL/ELE Obtained here by simulation for GDL thick- 
nesses from 150 to 350 um was in the range of 1-7 Q cm?, which 
is far too high and thus these results are unreliable. If the con- 
tact resistance really was that high, the voltage of a unit cell 
would drop to zero only due to resistive losses already at the 
current density of approximately 500 mA cm7? with the GDL 
compressed to 250 pm and Re GpL/ELe Of 2.5 Q cm?. A possi- 
ble explanation for these unrealistic results is that the current 
profile through the electrode and GDL interface was signifi- 
cantly uneven. Practically all of the current concentrated to the 
very edge part of GDL since the graphite and GDL conductiv- 
ities were significantly higher than the measured conductivity 
of electrode (320+ 20S m~'). This may have caused numer- 
ical problems in the simulations and distorted the results. In 
addition, due to the low conductivity of the electrode, the sep- 
aration between the graphite current collectors had to be kept 
very small, which led into inaccuracies in the assembly proce- 
dure of the measurement setup. However, the presented method 
should be applicable for the evaluation of the contact resistance 
between GDL and electrode, and more effort has to be put in 
order to achieve reliable values for it. 


4. Summary and conclusions 


The aim of this study was to evaluate the physical properties 
of GDL affected by compression in order to provide more insight 
into the mass and charge transport phenomena occurring in PEM 
fuel cells. The local properties of GDL may vary significantly 


due to inhomogeneous compression caused by the channel/rib 
structure of the flow-field plate. However, practically no previous 
experimental or modeling studies were found which properly 
consider inhomogeneous compression and its effects. 

It was observed that GDL under the channel remained almost 
at the initial thickness regardless of the width of the channel. 
Slight compression of GDL under the channel was observed 
only at the highest compressions. The change in thickness of 
GDL under the rib was most probably due to the loss of porosity, 
and hence in an assembled fuel cell the porosity of GDL varies a 
lot, causing significant local variation in mass, charge, and heat 
transfer properties. 

The in-plane gas permeability decreased non-linearly as the 
compressed thickness of GDL was decreased. For example, there 
is an order of magnitude difference in permeability between 
the uncompressed (380 um) and typically compressed (250 um) 
parts. The non-linear behavior may be due to changes in pore 
size and its distribution caused by compression. 

Compression had also a notable effect on GDL bulk con- 
ductivity. Both in-plane and through-plane conductivities were 
practically linear functions of compressed thickness. This was 
most probably due to the fact that the decrease in thickness led 
to loss of porosity, and the distances between conductive carbon 
fibers were decreased and more fibers came into contact with 
each other. 

As previously reported in the literature, contact resistance 
was found to vary significantly with the compression pressure. 
The difference in contact resistance between graphite and GDL 
was up to one order of magnitude between GDL compressions 
of 350 and 250 um. This result may be attributed to the increase 
in actual contact area of GDL at the interface as the compression 
pressure is increased. It is worth mentioning that the magnitude 
of bulk resistance of both GDL and graphite current collector was 
comparable to the contact resistance, on the contrary to the dis- 
cussions in previous studies. Thus, it is evident that the previous 
methods in which the bulk resistances were neglected and com- 
pression effects were disregarded may contain a considerable 
error in the prediction of contact resistance. 

Experimental observations in this work were made using only 
one type of GDL media, and thus quantitative results apply only 
on this material. However, due to the structural similarities, the 
trends detected should apply also on other carbon paper and even 
on carbon cloth GDLs. 

There are also other factors that affect the transport phenom- 
ena but are not evaluated in this study. These are the liquid water 
and heat transport properties. In an operating fuel cell, GDL 
operates in a wet environment because of the produced water. 
Liquid water transport may be one of the most decisive factors 
for the cathode performance as discussed, e.g. by Natarajan and 
Nguyen [44] or Lin et al. [45]. It is probable that the liquid water 
transport is significantly affected by the inhomogeneous com- 
pression and the consequent changes in pore size distribution. 
Phase change of water may be influenced by heat generated in 
fuel cell. Since interfaces could be one of the main heat sources 
in fuel cell as reported, e.g. by Noponen et al. [46], uneven 
distribution of electric and thermal contact resistances between 
GDL and electrode may create local hot spots. The hot-spots 
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may even lead into accelerated degradation of membrane, and 
thus shorten the life time of fuel cells. The evaluation of heat 
and water transport parameters is left for further research. 

Mass and charge transfer phenomena have a trade-off in terms 
of GDL compression. High compression increases GDL con- 
ductivity and decreases contact resistances, but hinders mass 
transfer. Thus, there may be significant differences in local 
cell performance. In order to gain more insight into the actual 
effects of the inhomogeneous compression of GDL, a modeling 
approach which takes into account these phenomena is needed. 
Such an approach is presented in the part II of this contribution 
[26]. 
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